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Abstract 

A relationship of the random walks on one-dimensional periodic lattice and the correla- 
tion functions of the XX Heisenberg spin chain is investigated. The operator averages 
taken over the ferromagnetic state play a role of generating functions of the number 
of paths made by the so-called "vicious" random walkers (the vicious walkers annihi- 
late each other provided they arrive at the same lattice site). It is shown that the 
two-point correlation function of spins, calculated over eigen-states of the XX magnet, 
can be interpreted as the generating function of paths made by a single walker in a 
medium characterized by a non-constant number of vicious neighbors. The answers are 
obtained for a number of paths made by the described walker from some fixed lattice 
site to another sufficiently remote one. Asymptotical estimates for the number of paths 
are provided in the limit, when the number of steps is increased. 
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1 Introduction 

The random walks is a classical problem both for combinatorics and statistical physics. 
The problem of enumeration of the paths made by the, so-called, vicious walkers on the 
one-dimensional lattice has been formulated and investigated in details by Fisher [1]. It 
is supposed that any two vicious walkers, provided both arrive at the same lattice site, 
annihilate not only one another but all other walkers as well. The problem mentioned still 
continues to attract considerable attention both of physicists and mathematicians [2]- [9]. 
Closely related problems arise also in the studies of the self-organized criticality [10], domain 
walls [11], and polymers [12]. In paper [13] a random walks of the annihilating particles on 
a ring was considered. In paper [14] a random turns walks on a semi-axes with a possible 
creation of the particles at the origin was studied. 

It has been shown in [15], [16] that the correlation functions, obtained as an averages 
over the ferromagnetic state of the XX Heisenberg chain, can be used for enumeration of 
the paths of random walks of vicious walkers. In the present paper the averages of special 
type are investigated both for the case of ferromagnetic state and for superposition of the 
eigen-states of the XX magnet in zero magnetic field. The averages in question play a role 
of the generating functions of number of paths of the vicious walkers. The calculation of 
the correlation functions is carried out by means of the functional integration [17], [18]. The 
answers are obtained for the number of paths of a single pedestrian which is travelling from 
one chosen site to another sufficiently remote lattice site. The asymptotical estimates are ob- 
tained for the number of paths in the limit, when the number of steps (and, correspondingly, 
the number of random turns) is increasing. 

The paper is organized as follows. Section [T] has an introductory character. The Hamil- 
tonian of the model and general calculation of the correlation functions are discussed in 
Section [2] Section [3] deals with the specific calculations and the corresponding asymptotic 
estimates. Discussion in Section [4] concludes the paper. 



2 The model and the correlation functions 

The XX magnet we are interested in is a particular limit of a more general spin model 
known as the XY Heisenberg chain, with the Hamiltonian in the transverse magnetic field 
h > given by [19], [20]: 

H = H + jH!-hS z , 

M 

n,m=l (1) 
. M M 

n,m=l n=l 

where S z is ^-component of the total spin operator, and 7 is the anisotropy parameter. The 
local spin operators cr* = (a^±ia^)/2 and a z n are given by the Pauli matrices, which depend 
on the lattice argument n G M. = {1, 2, . . . , M}, where M = (mod 2). The corresponding 
commutation relations have the form: 

*f] = 6 k>l af , [cr z k , of] = ±25 kil af . 
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The introduced hopping matrix is defined by the following entries: 



A nm = g ( 6 \n-m\,l + S<5|n-m|,A/-l) , (2) 

where 5 n j is the Kronecker symbol, and s can take two values: s — ±. It is assumed that 
the periodic boundary conditions o" +M = a™ are imposed for any n G M. The Hamiltonian 
if (ED) is reduced to the Hamiltonian of the XX magnet at zero value of the parameter 7. 

The most general definition of the time t and temperature T = 1//3 dependent correlation 
functions of the model under consideration looks as follows: 

G%{t) = ± Tr(^(0)of (t)e-? H ), Z = Tr(e^ H ), (3) 

where of (i) = e ltH a\e~ ltH and Tr means the averaging with respect to all eigen-states of the 
Hamiltonian H. In addition, the normalization involves the partition function Z. Calculation 
of the correlators (J3j) has been carried out in [21] as averaging over all eigen-functions of the 
Hamiltonian of the XX magnet. In [21] the main attention has been paid to a relationship 
between the correlation functions and the Fredholm determinants in the thermodynamic 
limit. In the present paper we shall consider the XX chain only and denote its Hamiltonian 
byH. 

To calculate the averages j3]) one can use a representation of the canonical Fermi vari- 
ables Cj, Cj, j G M. through the spin variables [19], [20]. The corresponding Jordan-Wigner 
transformation has the form: 



\ j=1 / v i=1 / 



where o| = 1 — 2cjcj. The periodic boundary conditions for the spin variables lead to the 
following boundary conditions for the Fermi variables: 

c M+ i = (-iyv 4 +1 = 4(-i)^ (5) 

where Af = ^2n=i c n c n is the operator of the total number of particles. The Hamiltonian 
H (PQ) takes the following form in the fermion representation [19], [20] 

H = H + P + + H-P-, (6) 

where P + (P~) are projectors on the states characterized by an even/odd number of fermions: 

P + + P_ = l, p + -p_ = (-if. 

The operators H^ 1 ([6|) are formally identical, their superscripts s — ± point out an appro- 
priate specification of the boundary conditions ([5|): 

CM+l = _sc li c l/+l = — s c \ ■ 

To put it differently, the quadratic in the fermion variables operators has the following 
representation: 

#± = c t£± c -^ 5 i/ ± = -Aw + / i /, (7) 
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where the matrices H ± are expressed through the hopping matrices (j2 
matrix: 

/ h -1/2 ±1/2 \ 

-1/2 h -1/2 

-1/2 h -1/2 



and / is the unit 



-1/2 



h -1/2 

\ ±1/2 -1/2 h J 

(only non-zero entries are displayed). Besides, the short-hand notations and c are used 
in ((71) for the M-dimensional row and column with the entries c^ n , c n , n ^ M.. 

In particular, the correlator at a — b — z takes the following form in the representa- 
tion (gD [18], [22], [23]: 



G%(t) = 1 



Trfcice 



f ttfoV**) + | Tr(c]c^ tH 4 Cl e-^ H ). (8) 



In order to calculate ((Hj), it is convenient to introduce the generating functional [18] 
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g = Q(S, T\\,v) = ^= Tr(e s e- XH e T e- uH ) 



(9) 



where A, v are the complex parameters, A + v — (3. The quadratic operators S = c^Sc and 
T = c^Tc, used in ([9]), are defined by means of the matrices S = diag{Si, S2, ■ ■ ■ , Sm} and 
T = diag{7i,T 2 , . . . ,T M }. For instance, the last term in right-hand side of (JHJ) is obtained 
from ([9]) in the following way: 

lim lim lim — — Q(S, T | A, u). (10) 
The trace in right-hand side of ([9|) can be re- written by means of ([6|) [18]: 



Ti{e s e- XH e T e- vH ) = \ (G+Z+ + Q~Z~ + Q+Z+ - G B Z B ), 



where 



and 



= Tr(e*e 

r?±<7± — rr^f ^-XH^ „T 



7 ± 

F 



Tr(e 



Tr(e e 



-l^e"^^ 



(11) 
(12) 



Tr((-1 



Moreover, for the partition function Z we obtain the representation: 

_ 1, _, _ . 



In the thermodynamic limit the terms with the subscript B are mutually compensated, 
therefore, in order to obtain Q j9]) it is enough to calculate Qp . 

The considered fermion representation is characterized by the existence of the Fock 
state |0) common for both operators H + and H~ , and satisfying the relations Ck\0) = 0, 
k e M. However, the corresponding coherent states over |0), 

M 



I z) = exp C clzk) I 0) = exp(c f 2:) | 0) 
v fe=i ' 

(z* \ = (0 1 exp 4 c k ) = (0 1 exp(2*c) 
v fc=i J 



are different for H + and H~ . Here the short-hand notations z* = {z\,...,z* M ) and z = 
(zi, . . . , zm) are used for the sets of independent Grassmann parameters Zk, z* kl k G M. (it 
is appropriate to omit the extra index ± in z*, z). Besides, Ylh=i c k z k = Elfeli ^ z k 
d::. Let us calculate GpZp in ( fT2l) using the representation of the trace in the Grassmann 
integration formalism [18]: 

g±Z± = j dzdz* e z * z (z*\e s e- XH± e T e- uH± \z). (13) 

In order to represent the right-hand side of this equality as the functional integral, let us 
introduce L new copies of the coherent states \x(I)), (x*(I)\, where I G {1,2, ...,L}. Each 
of the 2L multi-indices x*(I), x(I) is expressed by M independent Grassmann parameters. 
Using the decompositions of unity one can represent the right-hand side of ([13]) as the (L+l)- 
fold multiple integral. In order to express the quasi-periodicity condition it is convenient to 
introduce the auxiliary variables: 

-Ex{0) = x{L + l)=z, -x*(L + l)=x*(<d)E- 1 =z\ (14) 

where E = e s e~ XH± e T . Tending L to infinity, we obtain the functional integral over the 
space of the trajectories x*(r), x(r), where r e R: 

g±Z± = fe s d\*d\ Y[dx*(r)dx(T). (15) 

J T 

The action functional S = J L(r) dr is expressed through the Lagrangian L(r): 
L(t) = x*(t) (jt _ H^xir) + J*(t)x(t) + x*(t)J(t), 

where 

J*(r) = \*{6{t)I + 5{t - ^E- 1 ), J(r) = (S(r)i + 5{r - v)E)\. 

The integration over the auxiliary Grassmann variables A*, A in (fl~5l) guarantees the fulfilment 
of the constraints lfj"4"J ). The 5-functions in J*(r), J(r) reduce r G R to r G [0, 0\. The 
stationarity conditions SS/6x* = 0, 5S/5x = result in the following regularized answer [18]: 

5 - = det + J ■ (16 » 

Furthermore, we substitute (|T6l) into (llOp and pass to the momentum representation. The 
procedure described can also be applied to other correlators G^ t (t) Q, where a, b G {+, — }. 



3 Random walks 

As it has been shown in [15], [16], the flips of spins on a one-dimensional lattice may 
be associated with a random movements of walkers. Indeed, let us consider a state of 
the XX Heisenberg chain, which corresponds to the ferromagnetic ordering of M spins: 
Iff) = 0*Li |T)n (i-e., all spins are oriented "up"). Consider the average of the following type: 

^;/(A) = ^k+ e - AH Vrl^), (17) 
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where the notation Hq implies that the zero magnetic field h — is taken in the Hamilto- 
nian (jSJ), © (we shall omit the same subscript for the corresponding matrices H ± (J2J)), and 
A G C is an "evolution" parameter. "Up" (or "down") direction of spin corresponds to the 
empty (or filled) site. Differentiating Fj-i(X) (fTTj) and expanding the commutator [Ho, o\-~] 
we obtain the difference-differential equation: 

^F J ,(X) = l(F J+1 ,(X) + F^ 1 ,(X)) (18) 

(and similar equation can be also obtained for the fixed index j). Solution of the given 
equation is specified by the boundary conditions imposed on the lattice argument, and by 
the initial condition at A = 0. 

The average Fj-i(X) can be considered as the generating function of the trajectories with 
random turns that start at the Z-th site and end up at the j-th site. Indeed, let us introduce 
the notation T>% for the operator of differentiation of K-th order with respect to A at the 
point A = 0. The application of V% to the average (fTTjl leads to the answer: 

Vx [FjfiW] = (t \a+(-Ho) K ar\ ft> = £ . . . A^Ag. (19) 

ni,...,fix-i 

The right hand side of (TT9l) coincides with the entry at the crossing of the j-th row and the 
/-th column of the matrix given by the product of K copies of the hopping matrix ((2j) . Each 
matrix in this product corresponds to a transition between the two nearest sites of the lattice. 
After multiplication by 2 K (this is due to the accepted normalization of the matrix (j2J)), the 
right-hand side of (1T91) gives the number of the trajectories that consist of K steps and are 
connecting the /-th and j-th sites. Let us denote this number by \Pk(1 — > 

Let \Pk{W, ■ ■ ■ , In — > ji, ■ ■ ■ ,3n)\ be a number of trajectories consisting of K links made 
by iV vicious walkers in the random turns model. Here, the initial and final positions of 
the walkers on the sites are given respectively by the sequences h > h > ■ • • > In and 
ji > 3i > • • • > 3n- Let us consider the iV-point correlation function (N < M): 

F h j*,..^h,h,.,iM = (HWM ■ ■ ■ < e_Airo <V^ • • • ^lt>- (20) 

The present correlator is related to enumeration of the admissible trajectories which are 
traced by iV vicious walkers. Indeed, the application of the operator T>^ 2 to (1201) results in 
the average of the type 

(tk>S---<(-2^o)^;...%l^>- 

This average provides the numbers \Pnih-, ■ ■ ■ , In — * hi ■ ■ ■ ,3n)\ that can be established with 
the help of the commutator 

K 

[Ho, af^ . . . afj =5>7- ^ [Ho, <? h K k+1 ■ ■ ■ ° hi (21) 
fc=i 

(in this case, differentiation with respect to A/2, instead of A, allows to take into account the 
normalization of the hopping matrix (J2J)). The condition of non-intersection of trajectories 
of the walkers is expressed by the vanishing of the correlation function (l20l) for any pair of 
coinciding indices Ik or j p . 
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Differentiating ( !20l) with respect to A and applying pTT) . we obtain the equation: 



fe=l 

Equation (l22l) has been considered in [16] for the case of periodicity with respect to the 
lattice argument and with the initial condition: 

Fji,... ) j N ;h,...,l N \Q) 

The function Fj 1) j 2r __ t j m i l! i 2t _ j i N (X) can be expressed as the determinant of the matrix con- 
sisting of the averages of the type of (fTTj) [16]: 

F h ,..., jmh ,...,i N (X) = <!<•>(/•• : /jA)),. _ v . (23) 



N 




m=l 



3.1 Random walks on the axis 

Let us consider an infinite chain (M — > oo). Then, the modified Bessel function Ij-i(X) turns 
out to be a solution of equation (fl~8l) . which respects the condition -Fj ; /(0) = 8jj [15]: 

F i; ,(A) = Ij-i(X) = ±- r d0e AcoB V-^. (24) 

There exists the following expansion into the power series for 7j_/(A): 

1 /A \ Q 

where the summation index Q is subjected to the requirement: Q + \ j — l\ =0 (mod 2). In 
the limit of large "time" (A — >• oo) and for moderate values of m = \l — j\, using the known 
asymptotics for the Bessel function, we obtain for the generating function: 

n e A / 4m 2 - 1 \ 

where the decay is governed by the critical exponent £ = —1/2. 

Let the number K satisfies the relations K > — and K + \j — l\ =0 (mod 2). Then, 
differentiation of the series (1251) leads to the binomial relation \Pk(1 — > j)\ = Ck f° r the 
number of all lattice paths of the "length" K between two sites on the infinite axis: 

Here L denotes the one-half of the total number of turns: L = (K — m)/2. 

Let us consider now the multi-point correlation function Fj lj j 2i ...j N .i 1 j 2i ...j N (X). As it has 
been shown above, D\/ 2 [Fj u ...,j N ;i u ., .,i N W] has the sense of the number of trajectories of N 
vicious walkers each of which has made K steps. A different combinatorial interpretation 
of this object, however, can be proposed. Really, let us consider a representation of the 
multi-point correlator in the form of the determinant ( 1231 ). Its entries Fj r -i g (X) in the case of 
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an infinite chain are given by the Bessel function Ij r -i a (X) (|24l ). The operator T>f, 2 acts on 
the determinant as the differentiation of the product of N functions: 

(/i(x)/ 2 (x) • • • f N (x)) W = £ P(m, n 2 ,..., n JV )A (ni) / 2 (n2) • • • A niv) . (27) 

"l+"2H \-n N =K 

The notation J*™) = d n f[x)jdx n is used here, and the coefficients P(ni,n 2 , . . . ,n N ) are the 
numbers of permutations with repeats: 

P(m,7i2,...,7ij>r) = : — : : • 28 

Til! n 2 \ ... n^! 

Summation in (1271 ) is over all non-negative values of rii,n 2 , ■ ■ ■ ,tin, provided their sum is 
equal to K. 

Suppose further, that an iV-dimensional (hyper-)cubic lattice of infinite extension is given, 
and each site of this lattice is labelled by a set of N numbers. Let Tx(qi, q 2 , ■ ■ ■ , <7jv) be the 
number of the lattice trajectories that can be traced by some walker from the "initial" point 
O = (0, 0, . . . , 0) to a point (q%, q 2 , ■ ■ ■ , <7at) in K steps (by a single step the walker can move 
to one of the nearest sites). Let all numbers q^ be non-negative, and let the inequality 
K > qi + q 2 + • • • + qN be fulfilled, which means that the steps that can compensate each 
other are allowed. Let us denote the number of these steps as 2L, 



= K -q 1 -q 2 q N 

2 

Taking into account fl29l) . the following formula for the number of paths takes place: 



(29) 



T K (qi,q2, ■ ■ ■ ,gjv) = 2J P(qi + L 1} q 2 + L 2 , . . . ,q N + L>n,Li,L 2 , ■ ■ ■ ,L N ), (30) 

L 1 +L 2 +--+L N =L 

where summation is taken over all non-negative values of Li, L 2 , . . . , Ljv, provided that their 
sum is equal to L, and the formula (l28l) for the number of permutations with repeats is used. 

Turning back to the function Fj 1} j 2j _ } j m i 1 ,i 2 ,...,i N (\) let us define the matrix (n rs ) i< r ,s<N 
with the entries n rs = j r — l s . Then, we arrive to the following 

Proposition. The number of trajectories consisting of K links, which are traced by N vicious 
walkers on an axis, is expressed through the number of trajectories of the same "length" K , 
which are traced by a single walker travelling over sites of N -dimensional lattice of infinite 
extension: 

\P K (h, ■ ■ ■ , In — > ji, ■ ■ ■ , Jn) \ = ^ ) V2[^ii>-Jjv;«i,...,/jv(^)] = 

= ^2 (- 1 ) Vs TK{n ai i,n a22 ,...,n aNN ), ( 31 ) 

where summation is taken over all permutations S ai ,a 2 ,...,a N = S{a{, a 2 \ a N ) of the numbers 
1,2, . . . , N , and Vs implies a parity of a specific permutation. 

Proof. In order to verify (IBTT) one should develop the determinant (l23l ) by a row or by a 
column and then apply the induction using the relations (I26l) - (i30l . □ 
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Let us calculate, for instance, (1301 ) at N = 2: 

L 



T K (qi, q 2 ) - Cqltq2+2L C qi+L C q2+L ~ ^K^^Ki ( 32 ) 



k=0 



where L = (K — q\ — q 2 )/2 denotes one-half of the total number of turns. Then, using (I3TT) 
we obtain: 



v \/2[Fh,h-,h,h( X )] = 7k{nu,nz2) -T K (n 2 i,n 12 ) 



(33) 



where L = (K — n u — n 22 )/2 and the equality n n + n 22 = n 12 + n 21 is used. 

Representation of the entries l j 2) ...,j N ;i 1) i 2l ...,i N {\) (l23l) in the integral form (|24"1) allows 
to obtain the following expression [15]: 



*■* - -»(^)=wn(r§) e " AE " ,, ^ e * ,x 



(34) 



where S*(e m , e m , e idN ) is the Schur function [24], 

C ( \- det ( X 7 +N ~ k )l<J,k<N 

D- K {Xi,X 2 , . . . , Xn) = . N _ k , • \00) 

detixj )i<j, k <N 

The Schur function fl3"5l) depends on the partition 7r = (ttx, tt 2 , . . . , 7Ttv) defined by a sequence 
of non-negative integers, which are ordered according to non-strict decreasing: 7i"i > n 2 > 
■ ■ ■ > n N > 0. In virtue of translational invariance it is always possible to choose the numbers 
l\ > l 2 > ■ ■ ■ > In > —N for the initial position of the walkers and to define the elements 
of the partition by the equalities = h + k. In order to calculate the leading asymptotics 
of the generating function in the limit A — > oo, let us transform the integral (1341) into the 
following integral [7], [25]: 



/« n i%-^- A/!Ef - 9 >=^(n4 

l<j<k<N V=l ' 



It is a special case of the Mehta integral, which arises in the theory of the Gaussian matrix 
ensembles. Finally, we obtain the following asymptotics of the generating function for the 
trajectories traced by N vicious walkers: 



rlJ A^ 2 /2' ^ (2nW 2 H k-j 



where the well known formula for S\(l, 1, . . . , 1) is taken into account in A [9], [25]. There- 
fore, the power-like behavior of Fj x j N -i x i N (X) is characterized by the exponent £ = 
-N 2 /2, 
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3.2 Random walks over superposition of the eigen-states 

The eigen-functions of the XX Hamiltonian, given by the relations (EJ) , (EJ) , are constructed 
as combinations of the states, obtained by "flipping" of N spins in the state Iff) [21]. Indeed, 
let us consider all admissible strict partitions /x = (/ii,/i 2 , at), where M > ji\ > > 
■ ■ ■ > jiN > 1, and establish a correspondence between each partition and an appropriate 
sequence of zeros and unities: {e k = ek(i-i)} keM , where e k = <5fc lAtn , 1 < n < N . The required 
eigen-function is defined as: 

\y N { Ul ,...,u N )) = J2 T N({u k }\ »)(a M )^(a M _ i y^...(^r\ f), (36) 

{ek(fJ-)}keM 

where summation is taken over all strict partitions \x of the given type. The number of such 
partitions is expressed through the number of permutations with repeats (l28l ): P(N,M — 
N) = C^j. The wave functions satisfy the periodic boundary conditions, 

T N ({u k }\ii) = det^ 1 ) !< k ,i<N (37) 

are parametrized by the partitions fj, and by different, up to permutation, sets {u\, . . . ,un} 
of solutions of the Bethe equations: 

u\ M = (-lf-\ 1 < k < N. (38) 

These solutions have the form: u\ = e l27r/fe//M , where I k are integers or half-integers (this 
depends on parity of N). Due to the antisymmetry of (l36l) with the respect to permutations 
of the parameters u k , it is sufficient to restrict oneself to the strict partitions M > Ii > I 2 > 
• • • > Ijsr > 1 in order to guarantee the single- valuedness of u k . With the help of (l36l ) the 
corresponding normalized average 

« + ie- XHo ^) N ee -.. \"'"' (39) 

can be represented as a linear combination of (iV + l)-point generating functions f [20l) . There- 
fore, this average is related to the number of random walks of iV + 1 pedestrians. The initial 
and the final positions of one of them are fixed at h = 1 and j\ — m + 1, respectively, while 
for the rest (virtual) pedestrians the choice of their initial and the final positions is arbitrary. 

Calculation of equation (l39l) is of interest in the thermodynamic limit, when M and N 
are growing (their ratio remains finite), which means that the number of virtual pedestrians 
is increasing. In this limit [26] 

F m+1;1 (A) = « +1 e- XHo aY) N \ MN>>1 = TV'^+i^VT/), (40) 

where the notation Tr' points out that the procedure presented in Section [2] is used for the 
calculation of the normalized average. The difference-differential relation, analogous to the 
equation (JTHJ) is valid for F m+ i ; i(A) ( HOI ): 

^ F m+1;1 (A) = ~ (F m;1 (A) + F m+2;1 (A)) - Tr'(tf <T+ +ie - AH Vr) - 

-^'((^^)<e-^a^ (41) 





a m+l e 











10 



The form of the present equation makes it possible to suppose that the average F m+ i ; i(A) 
can also be of interest as a generating function of the random walks. 

Let us turn to the calculation of F m+1;1 (A) (14*01) in the fermionic representation (TJ)). It is 
convenient to reduce the problem to calculation of the generating function of the form: 

g = Tr'(e s c m+1 e- XH °c\e-» Ho ), (42) 

where the operator S is defined just like in (i.e., by means of the matrix S = 
diag{Si, S 2 , ■ ■ ■ , Sm})- Indeed, the functional -F m +i ; i(A) corresponds to the choice of v — 
and Sk = —iii9(m — k), where 6{m — k) is the Heavyside function, #(0) = 1. The second 
term in the right-hand side of ( 14T1 ) corresponds to the differentiation by v at the point v = 0, 
in the third term we put v = and differentiate with respect to S m+ i. In both cases we put 
Sk = —in9(m — k). Taking into account the fact, that the contribution of the terms labelled 
by the index B in (TTT1) is negligible at sufficiently large M and N, we approximately obtain: 



,u . d 



tr(e XH ei m+ i, 

da 



det(J + Mi + aM 2 )\ a=0 , 
M x + aM 2 = e-^°e § e- xS \i + agx^e"^), 



(43) 



where ei m+ i = (5i n 5 m+ i i)i< n) KM- The matrix H° is used instead of H ± since s can be 
replaced by zero for the sufficiently large M. 

The relation ([43]) is written in the coordinate representation. In order to pass to the 
momentum representation it is convenient to use certain formulas provided in [22]. Keeping 
the matrix notations as in (1431) . we obtain the answer for -F m +i ; i(A) (in the limit M — > oo, 
the corresponding operations should be understood in the sense of the operations over the 
corresponding integral operators [26]): 



F m+1;1 (A) = det(/ + W n 



tr(e- A£ "°e lim+1 ) - tr 



(44) 



(the notation tr, for instance, corresponds to the trace ofMxM matrices). The matrices 
U m , V m , ei jm+ i are given by the corresponding momentum representations of the matrices 
Aii, -M-2, ei, m +i ( |43l ). However, we shall need explicit expressions only for the traces tiU m 
and trV m (see below). In the momentum representation, i is a diagonal matrix of the 
eigen-energies of the XX model at h = [21]. Formally expanding F m+ i ; i(A) in the powers 
ofU m we shall obtain the answer in two lowest orders: 



-^m+i ; i(A) ~ F m+ i-i(X) + -F m +i ; i(A) trU m — tr V m , 

trU m = (M — 2m)F 1;1 (A), (45) 

m 

trV m = F m+ i ; i(2A) — 2 F m+ i ; /(A)F; ; i(A), 

i=i 

where the notation Fj-i(X) implies the relations ( l24l) . Although M and m are chosen to be 
large enough, the ratio m/M is assumed to be finite. Equation ( 1411 is fulfilled in each order 
separately by the terms presented in -F m+1;1 (1451) . 

By an analogy with the ferromagnetic case, let us act on F m+ i i(A) (1451) by the opera- 
tor T>^ 2 . Then, in the first order we shall obtain the relation (1261) . In the second order, the 
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answer is of the following form: 



MC L K C L K 



K m 



1=0 



1=1 



r L 



(46) 



By virtue of (l33l) . the result of the application of T>^ 2 to the second order function F^j^ ,/ 2 (A) 
is connected, as a particular case of (13TI ). with the number of the two-dimensional paths Tk, 
and is expressed through the corresponding determinant. It means that it will be appropriate 
to express (T4*6l) in the following equivalent form: 



(M - K)\P K {1 ^ I + m)\ 2 + 



1=1 



F m+ i-i(X) 



F m +l:l(X) 
Fl;l(X) 



K 

E 

1=0 



Fm+l-l (A) 

^ ; /(A) F t , L (X) 



(47) 



In other words, the result of application of Du 2 to (1451) in the second order can be refor- 
mulated in terms of the random walks of the two pedestrians (see ( 1231 ) and ( 1331) ) and the 
squared number of walks of a single pedestrian. The summation by the index I in tr V m (f4*5l) 
can be interpreted as the summation over positions of a virtual walker in (T4T1) . 

Using the equation (l33l) , one can represent (l47l) in terms of the number of trajectories on 
a two-dimensional lattice: 



771 

(M - 2(m + l))T K (m, 0) + ^ ^(m 



L L 



(48) 



In this relation various lattice trajectories of K links are enumerated. All these trajectories 
start at the same point O = (0,0) while they terminate on the segments of the dashed 
broken line which connects the points (L,L + m), (0,m), (m, 0), and (L + m,L) (see figure). 
Formally, the sign of the sum is not definite though its asymptotics is positive, in general. 
An analogous description is expected in the higher orders as well. 



m+L 



m 



O ' m m+L 

Typical configuration. 



Let us estimate the behavior of the number of paths, which is given by the representa- 
tion (T4*6l . in the limit, when the number of links K = m + 2L increases. We shall assume, 
that the restriction 1 < m C L is valid which means that m increases moderately in the 
comparison with the increase of the number of turns L: for instance, let L increase as m 2 . 
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Using the known asymptotical expansion of the logarithm of the gamma-function (see Ap- 
pendix) [27], one can estimate the binomial coefficient C\. Restricting ourselves by the first 
order of smallness, we obtain: 

" K „2/, A nf* m A m 2 \\ 2 K / m 2 \ ° 2L 



<?£ « * e -™ a /(*t) 1-— 1-— «-^= 1-— ~^=. (49) 



The second approximate equality in (1491 takes place if L is increasing faster than m 2 . The 
estimate (Jl9l) characterizes an increase of the number of the trajectories (l26l) for a single 
pedestrian. 

The third term in (l46l ) can be written as 2A(m, L)C^C^, where 

Mm, L) = -m + J2 1 (T , U ' • (50) 

Standard notation (a) n for the Pochhammer's symbol is used in ( |50l [27]. Applying again 
an expansion of the logarithm of the gamma-function (Al), we can estimate A(m, L) (1501 ): 

A(m, L) ~ mZ x {m, L) - Z (m, L) + C^ 1 ), (51) 

where 



2 

2 



m/2 
«=0 



Zi(m,I) = -l + e m2/4L -V 



m/2 



(52) 



m 

1=0 



Let the values m and L increase with the ratio L/m 2 being finite and of order of unity. It 
can be shown (by means of numerical check as well), that the coefficient functions Z (m,L) 
and Zi(m,L) remain finite in this case, and the contribution of Z (m,L) is negligible in 
comparison with mZi(m, L) in (I5T1) . One can use Eqs. (pf9l) and (I5T1) in order to estimate (1361) 
in the leading approximation: 

=— e- m2/2L (M + 2mZ 1 (m, L) - (irL) l/2 e m2/4L ). (53) 



Because of the behavior of the coefficient Zi(m,L), the corresponding contribution in (1531) 
may turn out to be comparable with M. The relation (l53l demonstrates that the description 
of the random walks considered in the representation of the superposition of the eigen-states 
is more complicated than the one in the ferromagnetic case. This description can be regarded 
as a simultaneous walks of the initial (i.e., principal) and virtual pedestrians. The ending 
points of the trajectories belonging to all the three segments of the dashed broken line 
on the figure (see the representation of two-dimensional random walks (1351) ) correspond 
to comparable contributions into the estimate (l53l) . In certain cases, characterized by the 
limiting behavior of the ratio m 2 /L, the contribution of the segment between the points 
(m, 0) and (0, m) can become dominating. 

4 Conclusion 

It is shown that the correlation functions of the XX Heisenberg magnet, calculated over the 
superposition of the eigen-states, as well as over the ferromagnetic state, are connected with 
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enumeration of the trajectories made by the walkers moving on the lattice. A relationship 
is established between the number of trajectories made by a several vicious walkers and the 
number of paths made by a single random turns walker on a lattice of a dimension equal 
to the number of the vicious walkers. Differentiation of the generating function, calculated 
over the superposition of the eigen-states, demonstrates a more complicated combinatorial 
picture than that of the ferromagnetic case. In particular, the set of the paths made by a 
single pedestrian is replaced by the set of trajectories made simultaneously by the principal 
and virtual (both vicious) pedestrians. An estimate is obtained for the number of trajectories 
made both by the principal and the virtual pedestrians. 
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Appendix 

Asymptotic expansion for the logarithm of the gamma- function at large \z\ and | argz| < n 
has the form [27]: 



\ogT(z + a) = (z + a - ^ log 2 - z + ^ log(27r) 



(Al) 



+ Et-ir 1 ^^ + o(— + 

^ v ' p(p+l) \z n + 

p=l v ' 

where n — 1, 2, 3, The Bernoulli polynomials B n (a) (Al) are defined as follows: 

n 

B n (a) = ^C^a""', 
1=0 

where C l n are the binomial coefficients, and Bi are the Bernoulli numbers. The first Bernoulli 
polynomials B p (a) look as follows: 

B (a) = 1, Bi(a) = a - ^ , B 2 (a) = a 2 - a + ^ , 

3 1 1 

B 3 (a) = a 3 a 2 + - a, B A (a) = a A - 2a 3 + a 2 , 

2 2 30 

where the Bernoulli numbers B\ are used: 

Bn = 1 , B\ — , Bo = — , Ba = . 

2 ' 6 30 
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